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1. Introduction 

In the last few years there has been an upsurge in interest in calculating the scalar matrix 
element on the lattice [jl]-||] either directly, by calculating the corresponding connected and dis- 
connected terms, or indirectly, via the Feynman-Hellman theorem. High statistics for dynamical 
simulations mean that reasonable signals can be obtained for disconnected terms and similarly the 
small statistical uncertainty on nucleon mass as a function of the quark masses enable reasonable 
fits to be made. Ideally, the results of both approaches should agree. 

Such calculations have also become particularly timely since the advent of the LHC because 
the scalar coupling fj = m q (N\qq\N) jm^ determines the fraction of the proton mass that is 
carried by quarks of flavour q. The strength of the coupling of the Standard Model (SM) Higgs 
boson or of any similar scalar particle to the proton is mainly determined by J^/r for q € {u,d,s}. 
Therefore, an accurate calculation of these quantities will help to increase the precision of SM 
phenomenology and to shed light on non-SM processes. 

The spin of the nucleon can be decomposed into a quark spin contribution A£ = Aw + Ad + 
As + ..., a quark angular momentum contribution L q and a gluonic contribution (spin and angular 
momentum) AG: 

l -= l -KL + L q + AG. (1.1) 

Experimentally, As is not well determined: HERMES obtained [|] As = -0.085(13)(8)(9) in the 
MS scheme. However, the signal is dominated by contributions in the small x region where models 
are used to extrapolate from the experimental results obtained at larger x. 

In these proceedings we present an update of an on-going project to calculate fj and Aq. In 
particular, the following improvements have been implemented since Lattice 2009 %: 

• Statistics on the 24 3 x 48 and 32 3 x 64 lattices have been significantly increased, by factors 
of roughly two and three, respectively, using the SFB/TR55 QPACE computers [§,|9|] (details 
are given in the next section). 

• An additional volume of 40 3 x 64 has been analyzed. 

The analysis is not yet finalized and we plan to further increase the statistics for the larger two 
volumes. The results presented here are preliminary. 

2. Simulation details 

The simulations were performed on rif = 2 configurations of nonperturbatively improved clover 
fermions with Wilson gauge action at j8 = 5.29 and fc sea = 0.13632. Details of the volumes and 
the number of trajectories analyzed are given in table ffl. The pseudoscalar mass corresponding to 
this fc sea value is around 270 MeV, using an inverse lattice spacing of 2.59 GeV determined from 
r (j8,K:) = 0.467 fm. 

The scalar and axial-vector matrix elements we are interested in are extracted on the lattice 
from the three-point functions corresponding to the diagrams given in fig. |l|. The axial-vector 
matrix element is related to Aq through, 

{N,s\qy^y 5 q\N,s) = 2m N sJ^-, (2.1) 
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Figure 1: The connected (left) and disconnected (right) diagrams associated with the scalar (r = 1) and 
axial vector matrix elements (r = JsJ,)- 



in Minkowski space notation, where is the nucleon mass and its spin (s 2 = — 1). The 
connected three-point functions were calculated separately by QCDSF, details of which can be 
found in refs. flip] 11]. The study presented here is only concerned with the disconnected terms. 



J3 = 5.29, n f = 2, fc sea = 0.13632 


Volume 


24 3 x 48 32 3 x 64 40 3 x 64 


# traj. 
m PS L 


^2000 ^ 1600 «730 
2.6 3.4 4.3 



For As and (A^|jj|7Y) only these terms contribute. 

We varied the quark mass of the current in- 
sertion (disconnected loop), as well as the mass of 
the valence quarks in the nucleon. In the follow- 
ing, we denote the K value corresponding to the 
quark loop and the valence quarks in the nucleon 
as Kj 00 p and fc va i, respectively. All combinations 
Kioo P ,Kvai G {0.13550,0.13609,0.13632} were used. 

These values correspond to the pseudoscalar masses, nips ~ 690, 440 and 270 MeV, respectively. 
The heaviest fc va i = 0.13550 roughly corresponds to the strange quark mass. 

The disconnected contributions to Aq and (N\qq\N) were extracted from the ratios of three- 
point functions to two-point functions (at zero momentum), 



Table 1: Details of the configurations used. 



Re(rgcg( ?f , ?i )£ x Tr(M- 



(x, t , x, f)ri OG p 



unpol 2pt 



ct{kA) 



(2.2) 



where the nucleon source and sink are at t\ and tf respectively, and the current is inserted at t. The 
three-point function is simply the combination of the nucleon two-point function, C2 p t(ff,?i), and 
the disconnected loop, £ x Tr [M^ 1 (x,?;x,?)ri oop ]. For the scalar matrix element we used, r2 pt = 
r un poi := (1 + 74)/2 and ri oop = 1. For Aq we calculated the difference between two polarizations: 
T2 P t = YjY5^unpo\ an d ri 00 p = 7; 75> wnere we average over all three possible j-orientations. 

In the limit of large times, ?f » ? 3> ?i, depending on the T-combination used, R dls will either 
approach the disconnected axial matrix element Ag dls or the disconnected scalar matrix element 
(N\qq\N) dis (once the vacuum contribution is subtracted). However, the statistical noise increases 
rapidly with increasing t—t\ and this time difference needs to be minimized, using smeared sources 
and sinks for the nucleon, in order to obtain a reasonable signal. A smearing study on a limited 
number of configurations indicated that the nucleon plateaued around t > 4a « 0.3 fm and we chose 
to insert the current at this timeslice (fj = 0). However, the higher statistics now available mean the 
ground state dominates around t >6a only. At zero momentum, the excited state contribution to 
R dls is governed by the time difference, tf — 1\, and we must be careful to choose tf large enough. 
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Figure 2: The results for {N\qq\N) dls and Aq dls extracted from the corresponding R dls as a function of the 
sink timeslice, ff, for all three volumes used and the heaviest fc va i = Kloop = 0.13550 combination. 

In fig. [2] we show the results for R dls as a function of tf > t for all three volumes studied and 
the heaviest KVal = Kloop = 0.13550 combination. So far we have chosen tf = 8a or 9a based on 
the quality of the plateau within given statistical errors. This depends on the observable and lattice 
volume. Once final statistics are reached we will fit the three- and two-point functions within R dl * 



of eq. ( |2.2D separately, as functions of tf, in order to extract the asymptotic values. 

The disconnected loop, £ x Tr [M _1 (x,r;x,?)ri oop ], was calculated using stochastic estimates, 
together with several noise reduction techniques: 



Partitioning [12, 13]: the stochastic source has support on eight timeslices. Additional two- 
point functions were generated for four time separated source points on each configuration. 
The forward and backward propagation from these four source points was combined with 
the loop to give us eight measurements of R dls per configuration. 



Hopping parameter expansion []14fl: for the clover action the first two terms in the expansion 
of the disconnected loop, 



Tr(M"T loop ) = 2icTr[(l - K^)- l T loop ] = Tr[(2Kl + 2k 1 I/> + K^M" 1 )!] 



loop J 



(2.3) 



vanish and hence only contribute to the noise. This means Tr [k 2 2 M~ l (x,t;x,t)ri oop ] can 
be used as an improved estimate of the loop. (In the case of ri oop = 1 the non- vanishing first 
term £ x 2fcTrll = 24kL 3 can easily be corrected for.) 



• Truncated solver method []15[]: 730 conjugate gradient solves were used, where the solver 
was truncated after 40 iterations. 50 BiCGStab solves running to full convergence were 
generated to correct for the truncation error. 

The noise reduction techniques other than time partitioning are only necessary for determining Ag; 
for the scalar matrix element the gauge noise dominates. 



3. The scalar matrix element: fj s and m q (N\qq\N) 



dis 



The results for fj s are presented in fig. || as functions of mp S corresponding to the mass of 
the (valence) quarks in the nucleon. No renormalization is required as the combination, m q (N\qq\N), 
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Figure 3: Results for the scalar matrix element: (left) fj s as a function of rrip S for the valence quark mass (i.e. 
the mass of the quarks in the nucleon) for the three volumes studied, (right) m q (N\qq\N) dls as a function of 
Wp S for the loop quark mass on the 32 3 x 64 volume, for the three valence quark masses. 



is scale and scheme independent. There is consistency between the values obtained on the different 
volumes for the heaviest valence quark mass, however, the spread between the results increases 
when the quark mass is reduced. Whether this is an indication of finite size effects will be clarified 
once the statistics for the 40 3 volume is increased. 

In fig. |5| we also display the results for the disconnected scalar matrix element, m q (N\qq\N), 
for the 32 3 volume as a function of mp s corresponding to the loop quark mass. This combination 
is relevant for extracting the sigma term, 

o N = m q (N\uu + dd\N) . (3.1) 

We found 2m q {N\uu\N) dls for fc va i = fq oop = fc sea = 0.13632 to amount to roughly 40% of the 
connected contribution to On for this volume. However, a more sophisticated method of extracting 
the matrix element from /? dls is required in order to make a firm comparison. 



4. The spin contribution: As and Aq 

The results for As on all three volumes are shown in fig. ||. No significant dependence on 
the valence quark mass nor on the lattice size is seen in the data. Neither is there any significant 
variation in the results if the loop quark mass is reduced. These numbers will have to be multiplied 
by a renormalization constant of approximately 0.8 for the MS scheme [16]. In contrast to the 
scalar case, the disconnected contributions are much smaller than the connected terms, at around 
10% for Ad and 5% for Au. 



5. Outlook 

In the short term, the aim is to reach our target statistics of 2000 trajectories for each volume. 
We then plan to begin an analysis close to the physical sea quark mass. The nonperturbative renor- 
malization for Aq needs to be calculated while for the scalar strangeness matrix element mixing 
with the light flavours needs to be considered. 
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Figure 4: Results for Aq: (left) As as a function of nip S for the valence quark mass (i.e. the mass of the 
quarks in the proton) for the three volumes studied, (right) Aq Als from the 32 3 x 64 volume for different loop 
quark masses, again as a function of nip S for the valence quark mass. 
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